Proximity effect in ferromagnet/superconductor hybrids: from diffusive to ballistic motion 
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We present an analytical study of the proximity effect in ferromagnet/superconductor (F/S) heterostructures, 
allowing for an arbitrary magnetic exchange energy as well as arbitrary impurity and spin-flip scattering rates 
within a quasiclassical approach. While previous studies mainly have focused on the clean or dirty limits, 
our results grant access to the regime of intermediate impurity concentrations, thus allowing us to probe the 
crossover from the clean to dirty limit. We find that in the crossover regime, all possible symmetry correlations 
of the proximity-induced anomalous Green's function are induced in the ferromagnet. We also point out that 
the local density of states oscillates spatially, not only for an F/S bilayer, but also for a normal/superconductor 
(N/S) bilayer in the diffusive limit, a fact which appears to have gone unnoticed in the literature. Within the 
weak-proximity effect regime, we present compact analytical expressions valid for arbitrary exchange fields 
and impurity scattering rates for i) the local density of states in an F/S bilayer, ii) the Josephson current in an 
S/F/S junction, and Hi) the critical temperature in an F/S/F multilayer. For all cases, we study in particular the 
crossover regime between diffusive and ballistic motion. Our results may be useful for analyzing experimental 
data in cases when the dirty limit is not fully reached, thus invalidating the use of the Usadel equation. 

PACS numbers: 74.25.Fy,74.45.+c,74.50.+r,74.62.-c 



I. INTRODUCTION 

The interest in ferromagnet/superconductor (F/S) het- 
erostructures has increased much during the last decad o 1 ' 2 ' 3 . 
This may probably be attributed to advances in experimen- 
tal fabrication/deposition techniques as well as intriguing the- 
oretical predictions. The main hope is that future devices 
and applications will rely on manipulation of not only the 
electron charge but also its spin. Based on this idea, a 
new research area known as superspintronics has emerged, 
aiming at utilization of charge and spin transport in ferro- 
magnet/superconductor heterostructures. For instance, sev- 
eral authors have investigated the possibility of dissipation- 
less currents of spin and charge in magnetically ordered 
superconductor s 4,5 ' 6 ' 7 ' 8 ' 9 ' 10 ' 1 A large number of other studies 
related to spin degrees of freedom in superconducting systems 
has also appeared in the literatur e) 12 ' 13 ' 14 - 15 . 

A considerable amount of attention has been devoted to the 
arguably most simple experimental laboratory where the in- 
terplay between ferromagnetism and superconductivity may 
be studied, namely a F/S bilayer. The two long-range or- 
der phenomena mix close to the interface, giving rise to in- 
teresting effects both from a basic physics perspective and in 
terms of potential applications. These effects include induc- 
tion of unusual superconducting symmetry correlations and a 
highly non-monotonic behaviour of various physical quanti- 
ties on the size of the system. The latter is a result of the 
non-uniform superconducting correlations that are induced in 
the ferromagnetic layer by means of the proximity effect. 

As a natural extension of the F/S bilayer, there has also been 
much focus on S/F/S systems and F/S/F systems, where the 
influence of ferromagnetism on the Josephson current and the 
critical temperature has been studied, respectively. The large 
majority of works related to these systems assumed that the 
diffusive limit was reached. In this case, elastic scattering 
on impurities renders the Green's function to be isotropic in 



space, while it may still retain a complicated spin-structure. 
From an experimental point of view, the diffusive regime is 
certainly relevant, but there are nevertheless some complica- 
tions. One point bears upon the theoretical framework used to 
study the physics in the diffusive regime. The quasiclassical 
Usadel 16 equation is widely employed to study the proximity 
effect in F/S heterostructures, and is valid under two main as- 
sumptions. Firstly, that the Fermi energy is much larger than 
any other energy scale and the essential physics is governed by 
fermions at Fermi level, and secondly, that the inverse impu- 
rity scattering rate is much larger than any other energy scale 
except for the Fermi energy. For strong ferromagnets such as 
Co or Ni, the second condition may be violated. In that case, 
one must revert to the more general Eilenberger— equation, 
which is only subject to the first condition. 

The Eilenberger equation is more complicated to solve ana- 
lytically than the Usadel equation, although some special lim- 
its permit fairly simple analytical expressions. Let h denote 
the exchange-energy of the ferromagnet while Tj mp denotes 
the inverse impurity scattering rate. The Usadel equation is 
then obtained from the Eilenberger equation by demanding 
h>Ti mp <C 1, while the case of a strong and clean ferromagnet 
is obtained in the limit hT\ mp 1. We assume that h ^> A 
is fulfilled. In Ref. 18 , some aspects of the DOS in F/S het- 
erostructures were considered to leading order in the param- 
eter (/iTi mp ) _1 , corresponding to a strong ferromagnet which 
falls outside the range of applicability of the Usadel equation. 
In Ref. 19 , the Josephson current in an S/F/S structure was also 
investigated for the case of a strong ferromagnet, hri mp ^> 1. 
Some authors have also considered F/S heterostructures where 
the impurity scattering rate was disregarded or assumed to be 
small, corresponding to the ballistic regime i 20 ' 21 ' 22,23,24 ' 25 ' 26 ' 27 

Although the agreement between theory and experiment in 
this research area has proven to be satisfactory in many cases, 
there are still discrepancies to be accounted for. For instance, 
the Usadel equation has failed to account quantitatively for 
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the critical temperature in F/S/F spin- valves. Furthermore, 
anomalous features in the DOS for a very thin F/S bilayer that 
could not be accounted for even qualitatively, were reported 
in Ref. 28 . Moreover, the Usadel equation approach fails from 
the start when addressing systems with strong ferromagnets. 

All of this points to the need of taking the role of impurity 
scattering more seriously. In this paper, we aim at doing pre- 
cisely so by solving the Eilenberger equation analytically and 
studying the crossover regime between ballistic and diffusive 
motion (see Fig. [TJ. To illustrate how various physical quan- 
tities behave in this crossover regime, we study i) the local 
density of states in an F/S bilayer, ii) the Josephson current in 
an S/F/S junction, and Hi) the critical temperature in an F/S/F 
multilayer for arbitrary values of h and r (within the quasi- 
classical approach). In each case, we present compact analyt- 
ical formula to facilitate comparison to experimental data in 
cases where the diffusive limit may not be fully warranted or 
where strong ferromagnets are involved. 
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FIG. 1 : (color online) Overview of the superconductor/ferromagnet 
heterostructure we will study in this paper. We take into account 
an arbitrary strength of the exchange field as well as an arbitrary 
rate of non-magnetic and magnetic scattering within a quasiclassical 
approach. 

This paper is organized as follows. In Sec. [TTJ we establish 
the theoretical framework which is employed in this work. In 
Sec. Hill we present our main results with belonging discus- 
sion: the DOS of an F/S bilayer in Sec. IIIIB1 the Josephson 
current in an S/F/S multilayer in Sec. IIII CI and finally the 
critical temperature in an F/S/F multilayer in Sec. HHP I Eqs. 
( [25) . (l32l . and (|37| | are the main analytical results of this work. 
We conclude in Sec. [IV] Below, we will use boldface notation 
for vectors, for 2x2 matrices, and . ". . for 4x4 matrices. 
The reader may consult the Appendix for a definition of the 
generalized Pauli-matrices we employ in this paper. 



II. THEORETICAL FORMULATION 

The Eilenberger equation reads 17 

iv F • V<? + [ep 3 + M — V imp - Sflip + A,g} = 0, (1) 

where g = g R (R, £,Pf) is the retarded part of the Green's 
function. Here, e is the quasiparticle energy, R is the center- 
of-mass coordinate, and pf (vf) is the Fermi momentum (ve- 
locity) vector. The self-energies that enter Eq. (UJ are the mag- 
netic exchange energy M = ft,diag{r3, T3}, the impurity scat- 
tering Vimp = — [ 1 /(2 T imp)](ff), the (uniaxial) spin-flip scat- 
tering 5 mp = — [i/(2Tflip)]/5 3 (g)p 3 , and the superconducting 
order parameter 

A=( 2 l ^ A ) 

All matrices used above {pt, n) are defined in the Appendix 
[Eq. (|4Tlll. The brackets (. . .) denote an angular average over 
the Fermi surface. Also, h is the exchange splitting while 
Timp(flip) is the scattering time associated with impurity (spin- 
flip) scattering. We may conveniently rewrite Eq. (UJ) as: 

iVf-Vfc + [(e + o-h)T3 + <tA + -i— (g^) 

+ ^-T 3 (g<r)T 3 ,g a ] = 0, a =T,|= ±1 (2) 

where the superconducting order parameter matrix A reads 

A= 0) ' A = A o eI *> ( 3 ) 

upon letting \ denote the phase corresponding to the glob- 
ally broken U(l) symmetry in the superconducting state. The 
brackets (...) denote angular averaging over the Fermi sur- 
face. We employ the Ricatti parametrization 29 of the Green's 
function: 

& = ^ C 7c ba 1 + u) • K = (1 + a ° b ° rK 

(4) 

Here, a a and b a are two unknown functions used to 
parametrize the Green's functions. They will be determined 
by solving the Eilenberger equation with appropriate bound- 
ary conditions. A general treatment of the Eilenberger equa- 
tion calls for a numerical solution. In the case of a weak prox- 
imity effect, however, the Eilenberger equation may be lin- 
earized in the anomalous part of the Green's function which 
permits an analytical approach. The assumption of a weak 
proximity effect corresponds mathematically to a scenario 
where higher order terms of {a a ,b a } are disregarded in the 
Eilenberger equation, i.e. one assumes that \a a \ <C 1, |6<t| <IC 
1. In an experimental situation, a weak proximity effect in 
F/S heterostructures may be expected whenever the tunneling 
limit is reached and the number of conducting channels at the 
interface is low. Also, assuming a superconducting reservoir, 
the proximity effect becomes weaker in magnitude upon in- 
creasing the thickness of the ferromagnetic layer. 



3 



The spatial depletion of the superconducting order param- 
eter near the S/F interface will be disregarded. This is an ex- 
cellent approximation in the corresponding low-transparency 
regime, which will be considered throughout this paper except 
for in Sec. MIDI where this issue is discussed further. At the 
S/F interface (x — 0) we use Zaitsev's boundary conditions. 
Define the symmetric and antisymmetric part of the Green's 
function as 

S<r,i = \{g<j,i + +ga : i~), A a ,j = \{g<y,i + -g<r,i~)> (5) 



where the ± superscript on the Green's function denotes 
right/left-going quasiparticle excitations and the subscript i 
denotes the ferromagnetic or superconducting region. The 
first of Zaitsev's boundary conditions 30 demands continuity 
of the antisymmetric part A a ,i of the Green's function. The 
second one relates the Green's functions in the ferromagnetic 
and superconducting regions to the interface transparency. We 
obtain 

A,,fW - A,, F 2 ) + j(5^s - S^f ) 2 ] 
T 

= -j[ Scr,F , Scr,s ]-> ( 6 ) 

where 1Z and T are the reflection and transmission coefficients 
satisfying 1Z + T = 1, and [...]_ denotes a commutator. High 
and low transparency interfaces correspond to T ~ 1 and 
T « 1, respectively. Although Eq. (0 is expressed rather 
compactly, a general solution for arbitrary T and 1Z is very 
hard to obtain. In the experimentally relevant situation, one 
may assume that T <C TZ. For a low-transparency barrier and 
a weak proximity effect, Eq. (|6]l simplifies greatly to 

Aa,F\x=0 — l[Sa,F,S a ,s]~\x=0, (7) 



where 7 = T/(41Z) is a measure of the barrier trans- 
parency. At the end of the ferromagnetic layer, we demand 

A a , F\x=d = Q. 

We consider here an effective one-dimensional calculation, 
which should provide sound results due to the isotropic na- 
ture of the ferromagnetic and superconducting order parame- 
ters. We do not expect any qualitative differences from a two- 
dimensional or three-dimensional model, since the supercon- 
ducting gap and the magnetic exchange field do not depend 
on the quasiparticle momenta, and since there are no surface- 
bound states 31 at the interfaces of the systems we consider. 
Thus, it should be possible to capture the essential physics by 
studying an effective one-dimensional model, which permits 
us to proceed analytically. This point of view is supported 
by the fact that, as seen later in this work, we reproduce in 
limiting cases previous results obtained in the literature which 
employed a two-dimensional calculation. 

Under the assumption of a weak proximity effect, the Eilen- 



berger equations in the ferromagnetic region take the form: 

aiv F d x a a + 2a a (e + ah) + -^—{a% - a~ a ) 

+ -^(3a« + a- Q ) = 

cav F d x b a - 2b a {e + ah) - ^—(b° - b~ a ) 

--^(3b« + b- a )=0, (8) 

^Tsf 

where a = ± denotes right- and left-going quasiparticles, re- 
spectively. It is necessary to take into account the direction of 
the quasiparticles at Fermi level due to the term vp • Vg in 
Eq. (Q~|i. Thus, a denotes the spin direction while a denotes 
the direction of motion in a% and likewise for b". The im- 
purity and spin-flip scattering self-energies enter Eq. (JHJ by 
means of the matrices V[mp and Sfl; p in Eq. ([1), which both 
depend on the Fermi-surface averaged Green's function. For 
a weak proximity effect, we have 

For a bulk ferromagnet, the solution is a^r = 6^ = 0. 

In Ref. 18 , the DOS in a S/F bilayer was studied by neglect- 
inb both spin-flip scattering (r S f — > 00) and the coupling term 
between the right- and left-going excitations in Eq. I©. In this 
case, one finds that Eq. <[8J reduces to 

±w F d x a± + [2(e + ah) + — ?— ]a± = 0, 

±w F d x bt - [2(e + ah) + — !— ]6± = 0. (10) 

The decaying solution for x — > 00 of the above equations 
reads 

at = K<7 exp[— K„x/l], b~ = k ba exp[— «x/Z], 

K CT = 1 - 2i(e + ah)r impi I = v F n m p. (11) 

while a~ = 6+ = 0. Above, k ac and are constants to 
be determined from the boundary condition at x = 0, and the 
structure of the Green's function becomes 

S„t - (I S) , A„, F - (_?_ < ) . (.2, 

This shows how the decay length of the proximity-induced 
anomalous Green's function in the ferromagnet is governed by 
the mean free path I, and that it is independent of the exchange 
field in this main approximation. We now present a more rig- 
orous solution by fully taking into account the coupling-term 
in Eqs. ®. To solve this problem, we note that Eqs. ([8j may 
be written as a matrix differential equation: 

d x a a — M aCT a CT , a CT = [a+a o :] T , 

M aa = — ( r ° 9 ), (13) 
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where T denotes matrix transpose and we have defined the 
auxiliary quantities 



r a = 2i(e + ah) - (g imp + 3g Bt )/2, 

? = (ffimp - 5sf)/2, ffimp(sf) = 7"imp( s f)- 



(14) 



Diagonalizing M a according to Da_ = Pa_ 1 M a P rT , we ob- 
tain the trivial set of decoupled differential equations 



9x&-a — DgJ^a; &l 

From the above, we find that 

5± = C± a e ±x ° x , A CT = ^V 
while the diagonalization matrix P a reads 



9 2 , 



(15) 



(16) 



p = {V\o P2a 
— \P2a Pla 



Pic 



Na: P2o 



, G a = g/(vFK + r a ), 
-N a G a , N a = (1 + |G ff | 2 )- 1/2 . 



(17) 



In the superconducting region, we employ the bulk solution 
under the assumption that the interface transparency is low 
and that the ferromagnetic layer is much more disordered than 
the superconductor 1 . In this main approximation, we may em- 
ploy the bulk solution of the Green's function in the supercon- 
ductor: 



± _ I c{6) as{9) 
// " " I -as{6) -c(9) 



with the definitions c{9) = cosh(6>), s(6) = sinh(6»), 6 = 
atanh(A/e). Once the expression for the Green's function in 
the ferromagnet has been obtained, one may calculate various 
physical quantities of interest. By approximating S a ,F — T3 
in Eq. (O in accordance with a weak proximity effect, we 
obtain for the case where the impurity-scattering coupling be- 
tween the Ricatti-equations is ignored: 

9%f = Ts+2jas(0)exp(-K rr x/l)(Ti ±172), (18) 

which is precisely the result of Ref. 18 for two semi-infinite 
superconducting and ferromagnetic layers in contact. When 
the coupling is properly taken into account, in addition to the 
vacuum boundary condition at x = d, we find that 



Cler = 
C 2 a = 



1 2a± 
26± -1 



upon defining 



piC la {X a )+p^C 2 a{Xa), 



2jas(6)e ) ' 



Pa - Pa 



2sinh(A CT d) 



jcrs(8)e 



\ a (d—x) 



(pt - Pa ) sinh(A CT d) ' 



(19) 



and p^ = Pi.2a- Note that in the diffusive limit where 
.9imp S> {h, e, Aq, g s f}, one would expect that the distinction 
between right-going and left-going particles is removed, such 
that g a .F + = 9a, f - This is easily shown by exploiting 



lim (v F \a + r a ) = -gimp/2, (20) 

p >{/>,e,A ,te} 



as seen from the previous equations. We also want to com- 
pare the results for g; mp ^> {h, e, Ao, <? s f} with those obtained 
when using the linearized Usadel equation. The Usadel equa- 
tion in a diffusive ferromagnet then reads 



Dd 2 J ± + 2i(e + ig st ±h)f ± = 0, 



(21) 



where f± — /t±/ s and f t is the odd-frequency triplet anoma- 
lous Green's function while f s is the even-frequency singlet 
anomalous Green's function (both are isotropic in momentum 
space). We obtain that the only physically acceptable (decay- 
ing for x — > 00) solution is 



f+=f e lk + x ife>0, /_=/ e- 
fc± = A /2i(e + i 5sf ± h)/D, 



ik— x 



if s < 0, 



(22) 



where fo is a constant to be determined from the boundary 
conditions. Above, D is the diffusion constant. For consis- 
tency, we should be able to obtain the same decaying solution 
from Eq. $1% when g lmp ^> {h, e, Ao, <? s f}- Focusing on the 
wavevector, we see that in this limit: 



-2i(e + ah)g imv + 2g imv g sf 



= yj-2i(e + ah + ig sf )/D, 



(23) 



where D = w^Ti mp is the diffusion constant in one dimen- 
sion (in three dimensions, D = ^Ti mp /3). Eq. (l23l is then 
consistent with the form of Eq. (l22l . 

With a complete description of the behaviour of the Green's 
function in the ferromagnetic region, we now investigate the 
influence of the proximity effect on the local density of states 
(LDOS), and also study the singlet and triplet superconduct- 
ing order parameters induced in the ferromagnet. The normal- 
ized LDOS as obtained from the solution of the Eilenberger 
equation may be written as 

N(e, x) = \ ^<Re{l + 4a CT (e, x)b a (e, x)}) (24) 



for a weak proximity effect. In the normal state, the normal- 
ized DOS is Nq = 1, Inserting the expressions for and b J 
into the above equation yields 
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N(x,e)=1--Rb\ 



2 7 V(0) 



smh(\ a d)(l + G a y 



1 + G\ - 2G CT 2 sinli(A CT d - 2\ a x) + 



r 



cosh(2A CT a;) 
sinh(A CT G?) 



(25) 



Eq. (l25l l is the first of our three main analytical results in this 
work. Within the weak proximity effect regime, it provides 
a general expression for the DOS, taking into account an ar- 
bitrary exchange field and impurity scattering rate. As seen, 
the correction to the normal state DOS Nq = 1 is zero for 
a vanishing interface transparency (7 = 0). While the weak 
proximity restriction only allows access to variations from the 
normal-state of DOS of around 10%, this seems to be suffi- 
cient for the experimentally relevant situation. For instance, 
the deviation from the normal-state DOS due to the supercon- 
ducting proximity effect was of order 1% in Ref.[32l 

In order to study the superconducting correlations inside 
the ferromagnetic region, first note that the full structure of 
the retarded Green's function is 



9 L 

-f -9 



where the spin-structure reads 



I = un fn 
- \fu fa 



(26) 



(27) 



and we have defined f a p — / q /3(pf, £> x) and /(pf, £, x) — 
/(— pf, — e, x)*. From the Ricatti-parametrization, we may 
define the different symmetry-components of the anomalous 
Green's functions as follows: 

/BSE = a ( a ° + <V )> /° S0 = a ( a i ~ O. 

/eto = £>+ - a-), /ote = £>+ + O- (28) 



Here, the abbreviations are explained in Tab. Q] Note that in 
the general case of finite h and Tj mp , all possible symmetry 
components of the anomalous Green's function are induced 
in the non-superconducting region. In the case of h = 0, 
one may confirm from Eq. ([T9T > that — > cra^, where 
is independent of a, such that /ote = /eto = 0. Physi- 
cally, the induction of other symmetry components than /ese, 
corresponding to the bulk superconductor, may be explained 
as follows 33 '-. In a normal metal/superconductor junction, 
the translational symmetry is broken at the interface separat- 
ing the two regions. This causes even-parity and odd-parity 
components of the Green's function to mix near the inter- 
face. Since the Pauli-principle must be satisfied at all times, a 
change in the parity-symmetry of the Green's function must 
be accompanied by a change in either spin- or frequency- 
symmetry. In the absence of an exchange field, nothing breaks 
the spin symmetry, such that only the frequency-symmetry 
may be altered indirectly by the broken translational symme- 
try. However, if the spin-symmetry is also broken by replac- 
ing the normal metal with a ferromagnet, the spin-symmetry 
of the Green's function may also be altered. These considera- 
tions are summarized in Tab. U The possibility of a bulk odd- 
frequency superconducting state was discussed in Refs . 35 ' 36 , 
and there has very recently been made some predictions con- 
cerning characteristic transport properties of such a bulk odd- 
frequency superconducting state . 33 ' 37 ' 38 ' 39 . 



TABLE I: Proximity-induced anomalous Green's functions in a normal metal in contact with a conventional BCS-superconductor, which has 
an even-frequency spin-singlet even-parity symmetry. Below, the quasiballistic limit regime is characterized by a vanishing or small value of 
the impurity scattering rate, while the diffusive limit is characterized by an impurity scattering rate which dominates all other energy scales in 
the problem (except for the Fermi energy). 



Symmetry 


h ^ 0, quasiballistic 


h — 0, quasiballistic 


h 0, diffusive 


h = 0, diffusive 


Even-frequency spin-singlet even-parity (ESE) 


V 


V 


V 


V 


Odd-frequency spin-singlet odd-parity (OSO) 


V 


V 






Even-frequency spin-triplet odd-parity (ETO) 










Odd-frequency spin-triplet even-parity (OTE) 


V 









m. RESULTS AND DISCUSSION 
A. Anomalous Green's functions 



quasiparticle transport, and hence dependence of the different 



The linearized Eilenberger equations allow us to study the 
direct crossover from the diffusive to the ballistic regime of 
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FIG. 2: (color online) Plot of the proximity-induced anomalous 
Green's functions in the middle of the ferromagnetic region (xjd = 
0.5) using h/A = 15 and e/A = 0.5. 



B. Density of states 

To demonstrate the applicability of Eq. d25t , we study 
in particular how the DOS depends on the crossover from 
the ballistic (gi mp = 0) to the diffusive limit (gi mp ^S> 
{h,e,A .g sf }). We will fix 7 = 0.05 and h/A = 15 to 
model a realistic experiment, corresponding to a weak fer- 
romagnetic alloy like Cui-^Ni^ or Pdi-^Ni^. The setup is 
shown in Fig. [3] It is well-known that the DOS oscillates in 
space upon penetration deeper into the ferromagnetic region 40 
due to the presence of an exchange field, a feature which is 
robust both in the clean and dirty limit. However, the energy- 
dependence of the DOS in the presence of an arbitrary im- 
purity concentration has not received much attention so far. 
This is because most works concerned themselves with the 
simplified Usadel equation (diffusive limit) or the Eilenberger 
equation in the absence of impurities (clean limit). 



Superconductor 
Ferromagnet 
Insulator 




symmetry components on the impurity scattering. In the ex- 
perimental situation, one usually probes the DOS at the F/I 
interface x — d, although it in principle is possible to ob- 
tain a spatially resolved DOS in the entire ferromagnetic re- 
gion by using local scanning tunneling microscopy (STM)- 
measurements. Let us first focus on x = d and consider the 
ballistic limit in which case simple and transparent analyti- 
cal expressions may be obtained from Eqs. ((28). In the case 
h 7^ 0, we obtain 



/ese = ^[-2 7 s(c?)]/sinh(A ff d), /oso = 0, 

/eto = 0, /ote = X)[-2<rys(0)]/ smh{X a d) (29) 



Note that for h = 0, A CT becomes independent of a, leading 
to /ote = 0. At first glance, this appears to be in contradic- 
tion to Tab. Usince the odd-parity components are absent even 
in the ballistic limit. However, evaluation of Eqs. ((28) for 
x 7^ d reveals that these components are in general induced, 
as they should be. It is remarkable that the odd-parity compo- 
nents vanish exactly right at the F/I interface. In the presence 
of a finite exchange field h ^ 0, however, the odd-frequency 
component /ote survives at x = d, and its influence on phys- 
ical quantities such as the DOS may be directly probed there. 
These results suggest that in order to investigate the influence 
of the odd-frequency superconducting correlations /eto and 
foso, one would have to measure the DOS at several positions 
in the ferromagnetic region and not only at the F/I interface. 
In Fig. |2] we plot the different symmetry components of the 
anomalous Green's function in the ferromagnet and their de- 
pendence on the impurity level. 



FIG. 3: (color online) Setup for our study of the density of states. 

In Ref. [ijl corrections to the normal-state DOS as induced 
by the proximity effect were calculated under the assumption 
that /ir; mp ^> 1. This case corresponds to a ferromagnet where 
the exchange field is considerably larger than the self-energy 
associated with the impurity scattering. This may describe ei- 
ther a strong ferromagnet (one must still demand /i < Ep) or 
a weak ferromagnet with weak impurity scattering. Neither 
of these cases are possible to treat with the Usadel equation. 
In the present work, however, we do not impose any restric- 
tions on the parameter /iTi mp , which allows us to study the full 
crossover regime. This may be important in order to obtain a 
larger degree of consistency between theory and experimental 
data in the case when the diffusive limit is not fully reached. 

In Fig. 2k), we study the energy-resolved DOS for an in- 
termediate range of impurity scattering. As a measure of the 
junction width, we use the superconducting coherence length 
in the clean limit £s = Vf/Aq. To isolate the role of the im- 
purity scattering, we fix the junction width at d/(_s = 0.5. For 
a superconductor with up = 10 5 m/s and Ao = 1 meV, this 
corresponds to d ~ 30 nm, which is experimentally relevant. 
As seen, the DOS exhibits a slightly oscillating behaviour as 
a function of energy when the impurity scattering rate g\ mp is 
comparable in magnitude to the superconducting gap. This 
effect becomes more obvious for wider junctions d/£s ^ 1> 
and is attributed to bound states appearing in the ferromag- 
netic film. We discuss this in more detail below. As g; mp 
increases, however, the DOS becomes featureless for subgap 
energies although one may still observe an alternating pos- 
itive and negative correction to the zero-energy DOS upon 
increasing g lmp . In Fig. I4J3), we plot the spatially-resolved 
DOS at e = for various rates of the impurity scattering, in- 
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eluding the case when hr im? ~ 1. As seen, the oscillations 
of the zero-energy DOS are reduced with increasing impurity 
scattering. We have also investigated the effect of spin-flip 
scattering for an intermediate value of the impurity concen- 
tration. The spin-flip scattering, here taken to be uniaxial, 
is pair-breaking and thus suppresses the proximity-effect in- 
duced by the superconductor. This aspect agrees withRef.lH 
which found that both the triplet and singlet components are 
suppressed with uniaxial and/or isotropic spin-flip scattering. 
For other types of magnetic scattering, such as planar spin-flip 
or spin-orbit scattering, the singlet and triplet components are 
affected very differently^. 

The oscillations of the DOS in S/F junctions are usually 
attributed to the oscillating decay of the Cooper pair wave- 



function in the ferromagnetic region. In an S/N junction, this 
decay is monotonous, and hence one would not expect to see 
any oscillations in the DOS. However, we underline that the 
impurity scattering plays an important role in this respect. In 
the ballistic case g — ► 0, the proximity of the superconductor 
induces Andreev-bound states with well-defined trajectories 
which propagate in the normal part of the system. The statis- 
tical distribution of all possible trajectories is peaked at given 
lengths, typically at trajectories corresponding to the first and 
second reflection processes at the interface. As a result, the 
DOS in a clean S/N junction acquires oscillations both as a 
function of energy and coordinate inside the normal region as 
seen in Fig. [5] upon averaging over all possible trajectories. 
This effect is known as Tomasch-oscillations 42 . 




However, there is another point which has appears to have 
been overlooked in the literature: namely that the spatial os- 
cillations of the DOS in a S/N junction at finite energies do not 
vanish in the diffusive limit. Hence, the oscillating DOS as a 
function of distance penetrated into the non-superconducting 
region is not a feature pertaining uniquely to F/S junctions, 
as have been implied in some works 4 ^. To see this, we plot 
the spatially-resolved DOS both for a F/S and N/S junction 
in Fig. |6]in the diffusive regime. The curves are obtained by 
using the framework of Ref. 41 , and thus correspond to a full 
numerical solution of the Usadel equation without restricting 
ourselves to the weak proximity effect regime. The oscilla- 
tions of the DOS in the N/S case may be understood by noting 



that the induced superconducting Green's function in the nor- 
mal region has a finite center-of-mass momentum q = 2s/vf- 
This is typically much smaller than the center-of-mass mo- 
mentum acquired in a ferromagnet, q = 2h/vp, which means 
that the corresponding oscillation length is much larger, but 
still present. 

Having stated this, it should be noted that the oscillating 
nature of the anomalous Green's function does not necessarily 
imply that the critical temperature dependence or the Joseph- 
son current in N/S multilayers is non-monotonuous, e.g. dis- 
playing 0-7r oscillations, since the energy-dependence of the 
Green's functions is integrated out when obtaining the criti- 
cal temperature or critical current. For an F/S junction, on the 
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FIG. 5: (color online) Plot of the a) energy-resolved DOS at x = d and b) spatially-resolved DOS at e/Ao = 0.5 for several values of the 
impurity scattering rate. Here, the exchange field is set to zero, corresponding to a normal metal, and d/£ = 5.0. 



other hand, the Cooper pair wave-function may retain its os- center-of-mass momentum depends on the exchange field h. 
dilating character even after the energy-integration since the 
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FIG. 6: (color online) Plot of the spatially resolved DOS for a diffusive N/S and F/S junction, respectively. In both cases, oscillations of the 
DOS are seen at finite energies. We have here fixed d/£ = 3.0 and r = 0.2, using the notation of Ref. 41 (here, £ = ^/D/Ao while r denotes 
the barrier transparency). 



C. Josephson current in Fig. [7] Denoting the phase at the left (right) superconductor 

as +X (~ X)> me total phase difference is given by (p = 2\- 

We now evaluate the Josephson current in an SFS junction 
for an arbitrary impurity concentration, with a setup as shown 
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FIG. 7: (color online) Setup for our study of the Josephson current. 

The current through the junction is evaluated by 

NpSoevF 



lj = 



deUnh(0e/2)Tr{(p 3 e F (g R - g A ))} 

(30) 



under the assumption of equilibrium distribution functions. 
Here, Sq is the effective area of the contact through which 
the current flows, while (3 = 1/T is inverse temperature. Ex- 
perimentally, one measures the current that flows through the 
junction, corresponding to the x-direction here. We employ 
the following boundary conditions: 



A a .F \x=0 — l[^cr,F,S^g]-\ x=0 , 
A a ,F\x=d = — r y[S a ,F,S^ 1 ^ >t ]-\ x= d, 



(31) 



and approximate S a ,p = t 3 as in the previous section, in ac- 
cordance with our assumption of a weak proximity effect. Af- 
ter some calculations, we arrive at the following expression 
for the Josephson current: 



I J = 47 NpSoevpI c sin ip, with the definition 

s 2 (6»)(l-G (T )tanh(/3e/2) 



In = 



G a ) sinh(A CT <i) 



(32) 



The reader is reminded of the definitions 



Ga =9/{V r a-r +r, 

r a = 2i(e + ah) - (g im? + 3g s{ )/2, 

9 = (5imp - 3sf)/2, #imp(sf) 



.2 



T; 



imp(sf) ' 



(33) 



Eq. d32l is the second of our three main analytical results in 
this work. It is probably the most compact way of expressing 
the Josephson current for arbitrary exchange fields and impu- 
rity scattering rates within the quasiclassical framework. It is 
thus suitable both for the case of a weak ferromagnet (such 
as the alloy Cui-^Ni^), and for strong ferromagnets (like Co 
or Fe) regardless of whether they are clean or dirty. In exper- 
iments performed with such strong ferromagnets, where the 
exchange field may be of order 100 meV (3> T c ), the Usadel 
equation is not valid at the same time as the clean limit may 
not be fully reached. In this case, one has to use an expression 
valid for the crossover regime, which emphasizes the impor- 
tance of Eq. Q2t . 

Below, we will study how impurity scattering affects both 
the width- and temperature-dependence of the critical cur- 
rent, as well as its belonging 0-7T phase diagram. Bergeret et 
al^ investigated this in the limiting cases of hr lmv <C 1 and 



/iTimp 3* 1. while the majority of studies so far considered ex- 
clusively the limiting case of diffusive motion. We here pay 
particular attention to the crossover between the ballistic and 
diffusive sector, which has not been investigated previously. 
To model inelastic scattering, we add a small imaginary num- 
ber to the quasiparticle energy, e — > s + 16 where 5 = 10~ 3 . 
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FIG. 8: (color online) Plot of the critical current as a function of 
junction width d. We have used T/T c = 0.2. 

In Fig. [8] we plot the width-dependence of the critical cur- 
rent for a temperature T/T c = 0.2. As seen, increasing im- 
purity scattering suppresses the magnitude of the current and 
also reduces the oscillation length l osc . The dependence of 
the latter on impurity scattering is shown explicitly in Fig. [9] 
Using the Usadel equation, it is predicted that the oscillation 
length of the critical current in the dirty limit should depend 
on the impurity scattering rate like y / 7zlh 



L imp ' ° -y/'Hmp (for a 
discussion of the characteristic decay and oscillation lengths 
in the clean and dirty limit, see Tab. I in Ref. 2 ). We obtain a 
good fit with this in Fig. |9]when g lmp ^> A. For values of (?; mp 
comparable to A, however, the oscillation length saturates at 
a finite value. In the ballistic limit, the oscillation length is 
known to depend on the exchange field like \/h. We have 
also confirmed this for several values of h when g lmv ~ A. 

Also, one notes from Fig. [8]that the decay length of the cur- 
rent increases with the concentration of impurities. It should 
be noted that the measure £ used as a length unit in this con- 
text is independent of the impurity scattering rate, since we 
are using £ = vf/Aq. This way, we ensure that the effects 
observed are really due to the increased impurity scattering. 
If we for instance had used the mean free path l m f p = VFTi mp 
as a measure for the junction width, the scale would have been 
different for each value of (?; mp in Fig. [8] We also underline 
that the dirty limit condition is that £/7 m f p ^> 1, while the size 
d of the sample may be either smaller or larger than £ as long 
as that condition is fulfilled. 

In Fig. [10] we pay particular attention to the case hr lmp = 1 
which is inaccessible in the Usadel framework. As seen, noth- 
ing qualitatively new shows up in the d-dependence or the T- 
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FIG. 9: (color online) Plot of the oscillation length of the critical 
current as a function of the impurity scattering strength g mp . We 
have used T/T c = 0.2. 



dependence of the critical current as compared to the diffu- 
sive limit, although the decay rate is considerably lower. We 
also investigate how the 0-7T phase diagram of the Josephson 
junction is affected by impurity scattering. This is most con- 
veniently plotted in the d-T plane. In Fig. QT| one observes 
several features. First, it is clear that the area occupied by the 
and 7r phases, respectively, diminishes with increasing g lmp , 
in agreement with the shortened oscillation length of Fig. [9] 
Secondly, it is seen that thermal 0-7r transitions are practically 
speaking impossible to observe for scattering rates satisfying 
<?imp < h. As the scattering rate is increased, however, the 
thermal transitions become possible when g; mp 3> h, or equiv- 
alently /iT; mp <C 1. In this regime, the Usadel equation is valid 
and we obtain consistency with previous results. At all scat- 
tering rates, the width-induced transitions are possible. 
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FIG. 10: (color online) Plot of the d-dependence and the In- 
dependence of the critical current for gj mp /Ao = 15, corresponding 

tO /iTimp = I- 



D. Critical temperature 

Finally, we investigate an F/S/F layers where the critical 
temperature of the superconductor is sensitive to the relative 
orientation of magnetization of the two F layers. This effect 
is usually dubbed to a spin-switch effect in the literature. Our 
setup is shown in Fig. [12] Tagirov 44 was the first to point 
out the interesting opportunity to "activate" superconductivity 
simply by means of switching the direction of the magnetiza- 
tion in one of the ferromagnetic layers. Since then, a num- 
ber of works have elaborated on the spin-switch effect both 
experimentall y 4 ^ 46 ' 47 and theoreticall y 40 ' 48 ' 49 ' 50 . In particular, 
a convincing numerical approach was developed in Refii 9 -. So 
far, however, almost all theoretical works focused on the dirty 
limit, in which the critical temperature may be conveniently 
calculated by using the Usadel equation in the Matsubara fre- 
quency representation. Although the obtained results compare 
well qualitatively with experimental data, an unsolved factor 
so far is the discrepancy of two orders in magnitude of the 



predicted effect. Recently, it was proposed and investigated 51 
if an asymmetry in the interface transparencies of the F/S/F 
junctions could be responsible for this, in effect one of the in- 
terfaces was much less transparent than the other. The authors 
of Ref. 51 concluded that this was not the case. At present, the 
single ferromagnet F/S/F devices to have been examined so far 
have used strong ferromagnets, which falls outside the range 
of applicability of the Usadel equation 45,4 ^ 7 .. In light of this, 
it would be interesting to go beyond the usual treatment with 
the Usadel equation and solve the more general Eilenberger 
equation to investigate the role of the impurity scattering. 

A general analytical solution for arbitrary proximity effect 
and barrier transparency is hardly achievable, as pointed out 
previously. Nevertheless, it is reasonable to expect that one 
may capture the essential physics in the weak proximity effect 
regime. In order to calculate the critical temperature for the P 
and AP alignment, we assume that the temperature is close 
to T c , which allows us to write the Green's function in the 
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FIG. 1 1 : (color online) Phase diagram in the d-T plane for the 0-n 
transitions of the critical current for several values of the impurity 
concentration. 



superconductor as follows: 




(34) 



since limA->o c (#) = 1- For the normal part of the Green's 
function matrix, this means that (1 — a a b a )/(l + a^bc) ~ 
1, while for the anomalous Green's function one thus has 
2a a /(l + a a b a ) ~ 2a a . The self-consistency equation for 
the superconducting gap reads in general^ 



A 



-Tr 



{( 



Pi 



- V2 V 



de(g 



A > 0, (35) 



where A is the attractive interaction and g K is the Keldysh part 
of the Green's function. For an equilibrium situation [g K = 
(g R — g A ) tanh(/3e/2)] in the weak-proximity effect regime 
with a temperature very close to T c , this reduces to 



de tanh ( 

\2T r 



Once the anomalous Green's functions {a^,b^} have been 
obtained, one may solve Eq. d37l i numerically to obtain T c in 
the P and AP configurations. Using boundary conditions ex- 
plained below, we solve for the anomalous Green's functions 
in both the ferromagnetic and superconducting regions and 
obtain the following equation determining the critical temper- 
ature: 



/ c \ r „ . r r ±p±2iex/v F n Ti a (l _ pZaieds /v F \ ^ -, 

1 — N F X / detanhf— )e _1 1 - cos(2ex/^) - V Re(^ ^"77 L ] = 0, (37) 



with the cut-off energy u>, a = ±, and finally Eq. (|37| | is the third of our three main analytical results in 

this work. It gives an expression for the critical temperature 
= e ±A " dp — G]f ft e TX <> dp , in an F/S/F junction for arbitrary exchange fields and impurity 

R ± = e ±X^ dF _ G Right eT A^ _ (38) 
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scattering rates within the framework of quasiclassical theory 
in the weak-proximity effect regime. 

In order to find {a^ , }, we must introduce proper bound- 
ary conditions at each of the interfaces in the setup (Fig. [T2"l >. 
The left ferromagnet is assumed to occupy the region —dp < 
x < 0, the superconductor is located at < x < ds, while 
the right ferromagnet occupies the space d$ < x < d$ + dp. 
Thus, the ferromagnetic layers are assumed to have the same 
thickness dp while the superconductor has thickness ds. Due 
to the complexity of the problem, we will assume rigid bound- 
ary conditions at the superconductor/ferromagnet interfaces, 
which amounts to continuity of the Green's function. Al- 
though the low transparency limit is probably more realistic, 
it is reasonable to expect qualitatively correct results in this 
approach. Moreover, since we already assume a temperature 
close to T c , the proximity effect would be almost completely 
absent if we in addition incorporated tunneling interfaces. In 
general, high transparency interfaces cause a depletion of the 
superconducting order parameter near the interface, which 
means that one should strictly speaking solve for the spatial 
depletion of the gap self-consistently. In our approach, we 
do not incorporate this depletion since we are aiming for an- 
alytical results. A full numerical approach would, however, 
doubtlessly improve the accuracy of the results presented be- 
low, but at the prize of losing the analytical information. 

At the ends of the ferromagnetic layers, we impose vacuum 
boundary conditions. In total, the boundary conditions then 
read: 

x = -dp : Jffi = 0, 
x = : sgf. = go_, 

j Right 

x = d s ■ ga_ = g a p , 
x = d s + dp : A^ 1S p = 0, 

(39) 

After straight-forward calculations, we obtain an expression 
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FIG. 12: (color online) Setup for our study of the critical tempera- 
ture. 

for {a^,b^} in the superconductor. A few comments with 
regard to the expression Eq. d37t are in order. Firstly, it 
should be noted that the expression for the critical tempera- 
ture in Eq. d37b depends on the position x in the superconduc- 
tor through the spatial dependence of the anomalous Green's 
function. This dependence is of course artificial and a result 
of the approximations we have made in the calculations; in 
a real experimental sample, T c is a property for the entire 
layer and does not depend on the position in the superconduc- 
tor. The reason for why we obtain an artificial x-dependence 
in the expression for the critical temperature is because we 
have neglected the spatial modification of the order parameter 
A in the layer. Employing a fully self-consistent calculation 
would remove the spatial dependence of T c in the gap equa- 
tion. However, for thin superconducting layers ds/S, <C 1, 
our approximation is expected to be good. A similar proce- 
dure has been used in several other works which calculated 
T c by means of the Usadel equation. In those works, it was 
assumed that the anomalous Green's function in the super- 
conductor varied very little as long as ds/£, <C 1 was satis- 
fied, and hence one could ignore the spatial dependence of the 
Green's function once it had been found. More precisely, T c 
was evaluated in the middle of the superconducting region. In 
our case, we will use the same approximation since our ap- 
proach is analytical in nature. The main contribution to the 
integral in Eq. d37l i comes from energies e < A, for which 
the terms including the coordinate x on the right hand side of 
the equation change very little as long as ds/£, <C 1. We will 
focus on the difference between the critical temperature in the 
P and AP alignments, defined as 



AT C = T C AP - T c p . (40) 



We will normalize all temperatures on T°, which is the bulk 
critical temperature of the superconductor in the absence of 
a proximity effect. As demanded by consistency, the critical 
temperature approaches when dp — > 0. We choose the 
cut-off frequency asw/ Aq = 30. 



With the analytical solution in hand, we now present a study 
of the critical temperature in the P and AP configuration, in- 
vestigating in particular the role of impurity scattering. First, 
we plot the critical temperature as a function of ferromagnetic 
layer thickness with a fixed superconducting layer thickness 
of ds/t; = 0.03 in Fig. [13] Using a superconductor with 
£ = 200 nm, this would correspond to a thickness ds = 6 
nm. To ensure the validity of our assumption that the anoma- 
lous Green's functions vary little with x throughout the su- 



perconducting layer, we plot the critical temperature both at 
x/ds — 0.50 (symbols) and x/ds — 0.01 (dashed lines). As 
seen, the difference is neglible. From Fig. QjJ one may infer 
that the critical temperature in the P configuration goes to zero 
much faster than in the AP configuration as a function of the 
ferromagnetic layer thickness dp. This supports the notion 
that the antiparallell configuration favors superconductivity in 
the middle layer. The effect of impurity scattering is seen to 
suppress the critical temperature, in general. 
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FIG. 13: (color online) Plot of the critical temperature in an F/S/F structure as a function of the ferromagnetic layer thickness (If for fixed 
ds/S, = 0.03. Note the different scale for cLf in the middle panel. The symbols denote the result for x/ds = 0.50 while the dashed lines 
correspond to x/ds ~ 0.01. 



One may understand intuitively why the antiparallell align- 
ment is favorable compared to the parallell alignment, since 
the average exchange field cancels in the former case. Qual- 
itatively, our results are consistent with the monotonic decay 
found for a high barrier transparency when using the Usadel 
equation 2 . However, a more realistic scenario would clearly 
be to invoke low barrier transparency boundary conditions at 
the S/F interfaces. Due to the complexity of the problem upon 
including an arbitrary amount of impurities, we have used per- 
fectly transparent interfaces here as a first approximation. It 
would nevertheless be quite interesting to extend this formal- 
ism to low transparency interfaces to investigate the role of 
impurity scattering under those circumstances. Especially, the 
role of g; mp with regard to the re-entrant behavior of T c would 
be worth investigating. Our analytical results may serve as a 
basis for extending this formalism to low transparency inter- 
faces in the case of an arbitrary value for /iTi mp , as opposed to 
/iT; mp in the Usadel regime. 

IV. SUMMARY 

We have investigated various aspects of the physics result- 
ing from the proximity effect in ferromagnet/superconductor 
(F/S) bilayers. In contrast to previous works, which were 
limited to either the clean or dirty limit, we have taken into 
account an arbitrary scattering rate for both non-magnetic 
and magnetic impurities. This has allowed us to access the 
crossover regime from the ballistic to diffusive regime of the 
proximity effect. We have derived analytical formula for i) the 
proximity-induced DOS of an F/S bilayer, ii) the Josephson 
current in an S/F/S junction, an d Hi) the critical temperature 



of an F/S/F structure. Our results are valid for an arbitrary 
ratio of the parameter hr lmv , and are thus applicable both to 
weak ferromagnetic alloys as well as permalloys in either the 
diffusive or clean limit. 
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APPENDIX 

The Pauli-matrices used in this paper are defined as 




1 F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Rev. Mod. Phys. 
77, 1321 (2005). 

2 A. I. Buzdin, Rev. Mod. Phys. 77, 935 (2005). 

3 Yu.A. Izyumov, Yu.N. Proshin, M.G. Khusainov, Physics- 
Uspekhi 45, 109 (2002) 



4 M. L. Kulic and I. M. Kulic, Phys. Rev. B 63, 104503 (2001). 

5 F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. Lett. 86, 
3140(2001). 

6 I. Eremin, F. S. Nogueira, and R.-J. Tarento, Phys. Rev. B 73, 
054507 (2006). 



14 



7 Y. Zhao and R. Shen, Phys. Rev. B 73, 21451 1 (2006). 

8 M. S. Gr0nsleth, J. Linder, J.-M. B0rven, and A. Sudb0, Phys. 
Rev. Lett. 97, 147002 (2006); J. Linder, M. S. Gr0nsleth, A. 
Sudb0, Phys. Rev. B 75, 024508 (2007). 

9 Y. Asano, Phys. Rev. B 74, 220501 (2006). 

10 P. M. Brydon, Boris Kastening, D. K. Morr, and D. Manske, Phys. 
Rev. B 77, 104504 (2008). 

11 T. Champel, T. Lofwander, and M. Eschrig, Phys. Rev. Lett. 100, 
077003 (2008). 

12 A. Brataas and Y. Tserkovnyak, Phys. Rev. Lett. 93, 087201 
(2004). 

13 D. Huertas-Hemando, Yu. V. Nazarov, and W. Belzig, Phys. Rev. 
Lett. 88, 047003 (2002). 

14 F. Giazotto and F. Taddei, Phys. Rev. B 77, 132501 (2008). 

15 J. Linder, T. Yokoyama, Y. Tanaka, and A. Sudb0, Phys. Rev. B 
78,014516(2008). 

16 K. Usadel, Phys. Rev. Lett. 25, 507 (1970). 

17 G. Eilenberger, Sov. Phys. JETP 214, 195 (1968). 

18 F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. B 65, 

134505 (2002). 

19 F. S. Bergeret, A. F. Volkov, and K. B. Efetov, Phys. Rev. B 64, 

134506 (2001). 

20 M. Zareyan, W. Belzig, and Yu. V. Nazarov, Phys. Rev. B 65, 
184505 (2002); M. Zareyan, W. Belzig, and Yu. V. Nazarov, Phys. 
Rev. Lett. 86, 308 (2001). 

21 J. Kopu, M. Eschrig, J. C. Cuevas, and M. Fogelstrom, Phys. Rev. 
B 69, 094501 (2004). 

22 K. Halterman, O. T. Vails, and P. H. Barsic, Phys. Rev. B 77, 
17451 1 (2008); K. Halterman, P. H. Barsic, and O. T. Vails, Phys. 
Rev. Lett. 99, 127002 (2007); K. Halterman and O. T. Vails, Phys. 
Rev. B 66, 224516 (2002). 

23 J. Linder and A. Sudb0, Phys. Rev. B 75, 134509 (2007). 

24 A. I. Buzdin, L. N. Bulaevskii and S. V. Panyukov, JETP Lett. 35, 
178 (1982) [Pis'ma Zh. Eksp. Teor. Fiz. 35, 147 (1982)]. 

25 M. Fogelstrom, Phys. Rev. B 62, 11812 (2000) 

26 N. M. Chtchelkatchev, W. Belzig, Yu.V Nazarov, and C. Bruder, 
JETP. Lett. 74, 323 (2001). 

27 I. V. Bobkova and Yu. S. Barash, Phys. Rev. B 71, 144510 (2005) 

28 P. SanGiorgio, S. Reymond, M. R. Beasley, J. H. Kwon, and K. 
Char, Phys. Rev. Lett. 100, 237002 (2008). 

29 N. Schop ohl and K. Maki, P hys. Rev. B 52, 490 (1995); N. 
Schopohl, |cond-mat/9804064| 

30 A. V. Zaitsev, Sov. Phys. JETP 59, 1163 (1984). 



38 



31 C.-R. Hu, Phys. Rev. Lett. 72, 1526 (1994); J. Yang and C.-R. Hu, 
Phys. Rev. B 50, 16766 (1994). 

32 T. Kontos, M. Aprili, J. Lesueur, and X. Grison, Phys. Rev. Lett. 
86, 304 (2001). 

33 Y. Tanaka et al. , Phys. Rev. Lett. 99, 037005 (2007). 

34 T. Yokoyama, Y. Tanaka, and A. A. Golubov, Phys. Rev. B 78, 
012508 (2008). 

35 V. L. Berezinskii, JETP Lett. 20, 287 (1974). 

36 E. Abrahams, A. Balatsky, D. J. Scalapino, and J. R. Schrieffer, 
Phys. Rev. B 52, 1271 (1995). 

37 J. Linder, T. Yokoyama, and A. Sudb0, Phys. Rev. B 77, 174507 
(2008). 

Ya. V. Fominov, JETP Lett. 86, 732 (2007). 

39 J. Linder, T. Yokoyama, Y. Tanaka, Y. Asano, and A. Sudb0, Phys. 
Rev. B 77, 174505 (2008). 

40 I. Baladie and A. Buzdin, Phys. Rev. B 64, 224514 (2001). 

41 J. Linder, T. Yokoyama, and A. Sudb0, Phys. Rev. B 77, 174514 
(2008). 

42 W. J. Tomasch, Phys. Rev. Lett. 15, 672 (1965). 

43 A. Buzdin, Phys. Rev. B 62, 11377 (2000). 

44 L. R. Tagirov, Physica C 307, 145 (1998); L. R. Tagirov, Phys. 
Rev. Lett. 83, 2058 (1999). 

45 A. Yu. Rusanov, S. Habraken, and J. Aarts, Phys. Rev. B 73, 
060505(R) (2006). 

46 I. C. Moraru, W. P. Pratt, Jr., and N. O. Birge, Phys. Rev. Lett. 96, 
037004 (2006). 

47 I. C. Moraru, W. P. Pratt, Jr., and N. O. Birge, Phys. Rev. B 74, 
220507(R) (2006). 

48 A. I. Buzdin, A. V. Vedyayev and N. V. Ryzhanova, Europhys. 
Lett. 48, 686 (1999). 

49 Ya. V. Fominov, N. M. Chtchelkatchev, and A. A. Golubov, Phys. 
Rev. B 66, 014507 (2002). 

50 T. Lofwander, T. Champel, and M. Eschrig, Phys. Rev. B 75, 
014512 (2007). 

51 P. Cadden-Zimansky, Ya. B. Bazaliy, L. M. Litvak, J. S. Jiang, 
J. Pearson, J. Y. Gu, Chun-Yeol You, M. R. Beasley, and S. D. 
Bader, Phys. Rev. B 77, 184501 (2008). 

52 We here notations and conventions close to the ones of J. P. 
Morten, M. Sc. thesis, Norwegian University of Science and Tech- 
nology, 2003, which in turn are similar to J. W. Serene and D. 
Rainer, Phys. Rep. 101, 221 (1983). 



